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The two-loop radiative photonic corrections to Bhabha scattering are computed in the leading 
order of the small electron mass expansion up to the nonlogarithmic term. After including the 
soft photon bremsstrahlung we obtain the infrared-finite result for the differential cross section, 
which can directly be applied to a precise luminosity determination of the present and future e'^e~ 
colliders. 
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Electron-positron Bhabha scattering plays a special 
role in particle phenomenology. It is crucial for extract- 
ing physics from experiments at electron-positron collid- 
ers since it provides a very efficient tool for luminosity 
determination. The small angle Bhabha scattering has 
been particularly effective as a luminosity monitor in the 
LEP and SLC energy range because its cross section is 
large and QED dominated Jj. At a future International 
Linear Collider the luminosity spectrum is not monochro- 
matic due to beam-beam effects. Therefore measuring 
the cross section of the small angle Bhabha scattering 
alone is not sufficient, and the angular distribution of 
the large angle Bhabha scattering has been suggested 
for disentangling the luminosity spectrum • The large 
angle Bhabha scattering is important also at colliders 
operating at a center of mass energy y/s of a few GeV, 
such as BABAR, BELLE, BEPC/BES, DA$NE, KEKB, 
PEP-II, and VEPP-2M^ where it is used to measure the 
integrated luminosity ■ Since the accuracy of the the- 
oretical evaluation of the Bhabha cross section directly 
affects the luminosity determination, remarkable efforts 
have been devoted to the study of the radiative correc- 
tions to this process (see Q for an extensive list of ref- 
erences). Pure QED contributions are particularly im- 
portant because they dominate the radiative corrections 
to the large angle scattering at intermediate energies 1- 
10 GeV and to the small angle scattering also at higher 
energies. The calculation of the QED radiative correc- 
tions to the Bhabha cross section is among the classical 
problems of perturbative quantum field theory with a 
long history. The first order corrections are well known 
(see and references therein). To match the impressive 
experimental accuracy the complete second order QED 
effects have to be included on the theoretical side. The 
evaluation of the two-loop virtual corrections constitutes 
the main problem of the second order analysis. The com- 
plete two-loop virtual corrections to the scattering ampli- 
tudes in the massless electron approximation have been 
computed in Ref. where dimensional regularization 
has been used for infrared divergences. However, this 
approximation is not sufficient since one has to keep a 
nonvanishing electron mass to make the result compatible 
with available Monte Carlo event generators . Recently 



an important class of the two-loop corrections, which in- 
clude at least one closed fermion loop, has been obtained 
for a finite electron mass 0] including the soft photon 
bremsstrahlung 0- A similar evaluation of the purely 
photonic two-loop corrections is a challenging problem 
at the limit of present computational techniques j^]- On 
the other hand in the energy range under consideration 
only the leading contribution in the small ratio m'^/s is 
of phenomenological relevance and should be retained in 
the theoretical estimates. In this approximation all the 
two-loop corrections enhanced by a power of the large 
logarithm \n{m1/s) are known so far for the small angle 
Q and large angle 0, 0| Bhabha scattering while the 
nonlogarithmic contribution is still missing. 

In this Letter we complete the calculation of the two- 
loop radiative corrections in the leading order of the small 
electron mass expansion. For this purpose we develop the 
method of infrared subtractions, which simplifies the cal- 
culation by fully exploiting the information on the gen- 
eral structure of infrared singularities in QED. 

The leading asymptotics of the virtual corrections can- 
not be obtained simply by putting rrig = because the 
electron mass regulates the collinear divergences. In ad- 
dition the virtual corrections are a subject of soft diver- 
gences, which can be regulated by giving the photon a 
small auxiliary mass A. The soft divergences are canceled 
out in the inclusive cross section when one adds the con- 
tribution of the soft photon bremsstrahlung • Here we 
should note that the collinear divergences in the mass- 
less approximation are also canceled in a cross section 
which is inclusive with respect to real photons collinear 
to the initial or final state fermions (IS.) . This means 
that if an angular cut on the collinear emission is suffi- 
ciently large, 6cut ^ i/m^/s, the inclusive cross section 
is insensitive to the electron mass and can in principle 
be computed with = by using dimensional reg- 
ularization of the infrared divergences for both virtual 
and real radiative corrections like it is done in the the- 
ory of QCD jets. However, as it has been mentioned 
above, all the available Monte Carlo event generators 
for Bhabha scattering with specific cuts on the photon 
bremsstrahlung dictated by the experimental setup em- 
ploy a nonzero electron mass as an infrared regulator. 
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which therefore has to be used also in the calculation of 
the virtual corrections. Thus we have to compute the 
two-loop virtual corrections to the four-fermion ampli- 
tude A^"^^ {rrie, A). The general problem of the calculation 
of the small mass asymptotics of the corrections including 
the power-suppressed terms can systematically be solved 
within the expansion by regions approach [l^ . To get the 
leading term in m^/s we develop the method applied first 
in Ref. to the analysis of the two-loop corrections to 
the fermion form factor in an Abelian gauge model with 
mass gap, which is briefly outlined below. The main idea 
is to construct an auxiliary amplitude (me , A) , which 
has the same structure of the infrared singularities but is 
simpler to evaluate. Then the difference A'^^^ — A'-^^ has 
a finite limit SA^^'' as me, A ^ 0. This quantity does 
not depend on the regularization scheme for A^'^'^ and 
A^^'' and can be evaluated in dimensional regularization 
in the limit of four space-time dimensions. In this way 
we obtain A^^'^inie, A) = A^^\me, X) +6A^^') +0{me, A). 
The singular dependence of the virtual corrections on in- 
frared regulators obeys evolution equations, which imply 
factorization of the infrared singularities [ig. One can 
use this property to construct the auxiliary amplitude 
^(^)(me,A). For example, the collinear divergences are 
known to factorize into the external line corrections . 
This means that the singular dependence of the correc- 
tions to the four-fermion amplitude on nie is the same 
as of the corrections to (the square of) the electromag- 
netic fermion form factor. The remaining singular depen- 
dence of the amplitude on A satisfies a linear differential 
equation and the corresponding soft divergences ex- 
ponentiate. A careful analysis shows that for pure pho- 
tonic corrections yl^^''(me,A) can be constructed of the 
two-loop corrections to the form factor and the products 
of the one-loop contributions. Note that the corrections 
have matrix structure in the chiral amplitude basis |l8l |. 
We have checked that in dimensional regularization the 
structure of the infrared divergences of the auxiliary am- 
plitude obtained in this way agrees with the one given 
m Refs. HE!- Thus in our method the infrared diver- 
gences, which induce the asymptotic dependence of the 
virtual corrections on the electron and photon masses, are 
absorbed into the auxiliary factorized amplitude while 
the technically most nontrivial calculation of the match- 
ing term SA^^^ is performed in the massless approxima- 
tion. Note that the method does not require a loop-by- 
loop subtracting of the infrared divergences since only a 
general information on the infrared structure of the total 



two- loop correction is necessary to construct ^^^^(me, A). 
Clearly, the method can be adopted to different ampli- 
tudes, mass spectra and number of loops. For the calcula- 
tion of the matching term 5^'^' beside the one-loop result 
one needs the two-loop corrections to the four-fermion 
amplitude and to the form factor in massless approxi- 
mation, which are available [1,|23|. For the calculation 
of A^^\me, X) one needs also the two-loop correction to 
the form factor for A <^ TOe <C s which can be found in 
Ref. li^l as a specific limit of the result for an arbitrary 
momentum transfer. We independently cross-checked it 
by integrating the dispersion relation with the spectral 
density computed in Ref. The closed fermion loop 

contribution, which is included into the analysis |2lll22l|. 
can be separated by using the result of Ref. . Let us 
now present our result. We define the perturbative ex- 
pansion for the Bhabha cross section in the fine structure 
constant a as follows: a = J2'^=o (f )" "^^^ leading 
order differential cross section reads 



doW 
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or 
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where x = {1 — cos6')/2 and 9 is the scattering angle. 
Virtual corrections taken separately are infrared diver- 
gent. To get a finite scheme independent result we in- 
clude the contribution of the soft photon brcmsstrahlung 
into the cross section. Thus the second order correc- 
tions can be represented as a sum of three terms cr^^^ = 



VV 



(2) (2) 

(Tsv + cTss , which correspond to the two-loop vir- 



tual correction including the interference of the one-loop 
corrections to the amplitude, one-loop virtual correction 
to the single soft photon emission, and the double soft 
photon emission, respectively. When the soft photon en- 
ergy cut is much less than me, the result for the two last 
terms in the above equation is known in analytical form 
and can be found e.g. in Ref. ^^I- The second order cor- 
rection can be decomposed according to the asymptotic 
dependence on the electron mass 
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where the coefficients S^^'' are independent on nif.. The re- 
sult for the logarithmically enhanced corrections is known 
(see ^3 and references therein). For the nonlogarithmic 
term we obtain 



?(2) 



8Ci 



1 - x + x'' 



4 8 2 10 ^ 

X - X -\ x^ 

3 3 3 



x'' TT^ (-12 + 16x - 18x2 -I- 6x^) Mx) 



+ {2x + 2x^) ln(l -x) + (-3x + + 3x^ - 42;*) ln{x) + (-8 + 16x - Ux'^ + Ax^) ln(x) ln(l - x) 



3 



+ (4 - lOx + Ux^ ~ IQx^ + 4x^) \n{l - x) + {l - x + x^) (16 + 8Li2(a;) - 8Li2(l - x)) 



[1- X + X'' 



; V Jo •'0 y ^ 1^3 24 8 24 6 / Vl8 360 36 180 



109 



+ 1440"^ +1^'^"^'^ ^^"^ + 



93 231 279 , 93 , 

1 x X H x'^ 

8 16 16 16 



3 13 

2 + T"-12" 



2 11 J 



X n 



+ (12 - 12a; + 8a;2 - a;^) C(3) 



/67 5 2 



ln(a;) 



9 43 



a;-' In-' (a;) + 



+ (-6 + 6x - a;2 - 4a;3) C(3) + 



48^"^ 96^ 
21 45 



17 5 29 . 9 4 

—x^ H a;'^ x^ 



x'^ 43 



a; a;^ 5 ^ 19 

- H 1 x'^ H x" 1 tt" 

4 2 24 48 



ln2(x) 



12 96 



+ 7^^ In (^) + i 3a; + 3X-' + -a; - —x^ + —x^ tt 



7 



73 



15 



—X -x" + a;^ + 1 



6 12 



35 



15 



4 

7 25 , 13 4 

--X x^ -"^ 

6 24 



12 



24 



6 24 



ln(a;) + 6 - llx 



7 . 45 2 . 3 7 4 / 17 
--6a;+-a;^-6x3 + -a;^+ 



X ln^(l - x) + 
29 



a; I TT + 
1 



48 
3 

—a 

8 6 



-3 + ^x ~ ^-x^ + ^a;^ ) ln(x) + 
4 4 8 ' 



7 41 



-X + + 



9 3 Q 

1- + -a; + -a; 



, 29 49 , 5 , 74,,,, 

-4 H a; x^ + -x'' + -x* ln(x) 

6 12 6 8 ' ^ ' 



In-'(l-x) 



x^ ln2(x) 



71 



13 
16" 

1 3 

X H x~ 

32 4 48 



-24^ 



+ :^xM In'' (1 - x) + 1 8 - 16x + 24x2 _ j^g^s ^ g^4 ^ M _ 3^ _^ 3 2 ^ 5^ 
32 / ^ ' [ V3 4 6 



11 



4x'^ + x"* + 2 - —X + -x^ + — - X** ln(x) 



ln(x) 



s 24 

X X I TT" 

3 

15 



-6 + ^x 



4 + 9x - y x^ + 2x3 j 



25 



x^ - 5x3 ^ 23;-* ln(l - x 



ln(l - x) + (2 - 4x + 6x2 _ 4^3 ^ 22;*) Li2(x)|Li2 (x) + | - 8 



+16x - 24x2 + 16x3 - 8x^ 



2 4 x2 5 o 2 4 

3 + r+y-r +r 



6 - 8x + 9x^ - 3x3 



2 2 2 



+2x* ln(x) 



ln(x) 



x2 x3 



/ 31 x3 \ 

+ (10 - 14x + 9x2) + f -8 + llx - — x2 + — + j ln(l - x) 

X ln(l - x) + (-4 + 8x - 12x2 + 8x3 _ 4^4^ l^i^ix) + (2 - 4x + 6x2 _ 4^3 _^ 2a;'') Li2(l - x)|Li2 (1 - x) + 

3 

-5x2 2a;3 + (-4 - X + x2 + 2x3 - 23,4^ j^^^,) + (6 - 6x + x2 + 4x3) j^j^^ _ _^ | _ i|_ + 

+3x2 _ 5^3^ + (6 - lOx + 10x3 - 6x'') ln(l - x)] Li3 (1 - x) + ^-2 + ya; - yx3 + 2x^^ (, 



2 

6 + 5x 



7x - -x2 - 4x3 + 6x*j Li4 (1 - x) + \ -Q + 4x + -x2 - 7x3 j \^Y~~^ 



(3) 



where ((3) = 1.202057... is the value of the Rie- 
mann's zeta- function, Li„(z) is the polylogarithm, — 
[1 - In (x/ (1 - x))] In {scut/e), e = ^fsjl, and Scut is the 
energy cut on each emitted soft photon. In Eq.© /q""* 
stands for the n-loop nonlogarithmic coefficient in the se- 
ries for the form factor at Euclidean momentum transfer 

/« = ^2/12-1, 

(2) 15 ,43 2 59 4 7r2ln2 9 
^° = y+96" " 1440" "^"4^^^)-^^) 



As it follows from the generalized eikonal representation 
[2^. in the limit of small scattering angles the two- loop 
corrections to the cross section are completely deter- 
mined by the corrections to the electron and positron 
form factors in the t-channel amplitude. In the limit 
X — > Eq. (O agrees with the asymptotic small angle ex- 
pression given in |9Ll24|. which is a quite nontrivial check 
of our result. Note that Eq. ^ is not valid for very small 
scattering angles corresponding to x<m1/s and for al- 
most backward scattering corresponding to 1 — x<m'^/ 
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FIG. 1: Logarithmically enhanced (dashed line) and nonloga- 
rithmic (solid line) second order corrections to the differential 
cross section of the small angle Bhabha scattering as functions 
of the scattering angle for -^s — 100 GeV and ln(£cut/e) = 0, 
in units of 10~^. 

where the power-suppressed terms become important. 

The second order corrections to the differential cross 
section (a/7r)^dCT*^^Vd<7^°'' are plotted as functions of the 
scattering angle for the small angle Bhabha scattering at 
y/s — 100 GeV on Fig. |^ and for the large angle Bhabha 
scattering at y/s = 1 GeV on Fig. We separate the 
logarithmically enhanced corrections given by the first 
two terms of Eq. Q and nonlogarithmic contribution 
given by the last term of this equation. All the terms in- 
volving a power of the logarithm ln(ec«t/£) are excluded 
from the numerical estimates because the corresponding 
contribution critically depends on the event selection al- 
gorithm and cannot be unambiguously estimated with- 
out imposing specific cuts on the photon bremsstrahlung. 
We observe that for scattering angles 9<18° and 0^166° 
the nonlogarithmic contribution exceeds 0.05%, and in 
the low energy case exceeds the logarithmically enhanced 
contribution for 9<8°. 

To conclude, we have derived the two-loop radiative 
photonic corrections to Bhabha scattering in the leading 
order of the small electron mass expansion up to non- 
logarithmic term. The nonlogarithmic contribution has 
been found numerically important for the practically in- 
teresting range of scattering angles. Together with the 
result of Refs. 0, for the corrections with the closed 
fermion loop insertions our result gives a complete ex- 
pression for the two- loop virtual corrections. It should 
be incorporated into the Monte Carlo event generators 
to reduce the theoretical error in the luminosity determi- 
nation at present and future electron-positron colliders 
below 0.1%. 
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FIG. 2: The same as Fig. Q but for the large angle Bhabha 
scattering and -^s = 1 GeV. 
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